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Two-photon production in pp scattering is proposed as a means of studying virtual-pion emission. 
Such a process is complementary to real-pion emission in pp scattering. The virtual-pion signal is 
embedded in a background of double-photon bremsstrahlung. We have developed a model to describe 
this background process and show that in certain parts of phase space the virtual-pion signal gives 
significant contribution. In addition, through interference with the two-photon bremsstrahlung 
background, one can determine the relative phase of the virtual-pion process. 

PACS numbers: 13.30.-a, 13.40.-f, 13.40.Gp, 13.40.Hq, 13.60.Fz 



I. INTRODUCTION 

Near-threshold pion production in proton-proton scat- 
tering has a long history [|, [^]. More recently it has 
attracted much attention after precise data have become 
available from experiments at IUCF || and These 
showed that this relatively simple reaction is apparently 
poorly understood. Earlier works showed large discrep- 
ancy between the calculations and the data [|[ g[. Later 
several different mechanisms such as heavy meson ex- 
change 0, D , off-shell structure of the T-matrix |9|, [l0| , 
heavy meson exchange currents |TT[ , and approximations 
used in the calculation of the loop contributions Jl^| 
have been proposed to explain this problem in the ear- 
lier calculations. Calculations in a relativistic one- 
boson- exchange model |l3| and non-relativistic potential 
model 14 1 on the other hand appear to reproduce the 
data rather well. 

In this work we propose to extend the available kine- 
matical regime for neutral-pion production by investi- 
gating the process of virtual 7r° emission which can be 
observed through its two-photon decay. The interest in 
this process is manyfold. For example, the importance 
of off-shell form factors or the off-shell structure of the 
T-matrix [To) can be investigated when an extended 
kincmatical regime is available for measuring this reac- 
tion. In addition, studying virtual-pion production below 
the threshold for real-pion production in proton-proton 
scattering implies an important simplification in the de- 
scription of the process since the inelasticities, always 
present in the pion-nucleon scattering, are absent for vir- 
tual pions. More importantly, interference of two-photon 
production via a virtual pion with the background due 
to two-photon bremsstrahlung will determine the rela- 
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tive sign of the matrix elements. This will allow for a 
better insight in the underlying pion-production process. 
The sign is relevant regarding the discussion of the con- 
structive v.s. destructive interference of the higher-order 
diagrams in a field-theoretical approach . 

To describe the "background", two-photon 
bremsstrahlung, cross section we have developed a 
Soft-Photon Model (SPM) for two-photon emission. 

In this context SPM implies a covariant model sat- 
isfying gauge invariance, which obeys the proper low- 
energy theorem for small photon momenta (for example, 
for the one-photon bremsstrahlung the leading two orders 
in powers of the photon energy satisfy model-independent 
constraints Q). As a first step towards such an SPM, 
we develop in section II a new SPM for the single-photon 
bremsstrahlung amplitude. This novel SPM, based on a 
power-series expansion of the T-matrix, combines ideas 
of the two SPM's which are frequently used for single- 
photon bremsstrahlung in pp scattering: the original 
SPM [[t9[ which is directly inspired by the deriva- 
tion of the low-energy theorem for bremsstrahlung by 
Low jl8| , and a later one proposed in ref. , which has 
been very successful in reproducing the observed cross 
sections [ p2| , p3| . The important distinction of the new 
SPM from the existing two is that no explicit contact 
terms, or the so-called internal contributions, need to 
be introduced. This feature makes it the most suit- 
able model for developing the Two-Soft-Photon Model 
(2SPM) as discussed at the end of section II. Such a 
2SPM may also be used to calculate the background two- 
photon signal in the search for di-baryon states 0]. 

In this work, where the emphasis is placed on the fea- 
sibility of detecting the virtual-7r° signal, we have also 
employed a relatively simple covariant model to describe 
the pion-emission process. This model is discussed in 
detail in section III. The predictions of this model are 
shown to reproduce data on real-pion emission. 

In section IV explicit calculations are presented 
for two-photon production where both mechanisms, 
bremsstrahlung and virtual-7r° emission, are taken into 
account. The parts of phase space are indicated where 
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the second mechanism is relatively large. It is also 
demonstrated that interference between the two pro- 
cesses is very important. 



II. THE SOFT-PHOTON MODEL 

A starting point in an SPM description of 
bremsstrahlung in pp scattering is that the domi- 
nant -pole- contribution to the amplitude is derived from 
the Feynman diagrams where the photon is radiated 
off the external legs Jl^]. To this leading order, several 
higher-order, non-pole, terms need to be added which 
may correspond to meson-exchange, form-factor, and 
rescattering contributions. The observation made by 
Low, which is the essence of the low-energy theorem [[l8| , 
is that in any description for the amplitude which has 
the correct pole structure and is gauge invariant, in 
a power expansion of the amplitude, the leading two 
powers are model independently given by an expression 
involving only on-shell (i.e. measurable) quantities, such 
as the non-radiative NN T-matrix, charge and magnetic 
moment of the nucleon. 

In the formulation of a SPM description for 
bremsstrahlung this model independence of the leading 
contributions to the amplitude is exploited. In princi- 
ple, the T-matrix entering in each of the pole diagrams 
needs to be evaluated at different off-shell kinematics. 
In an SPM one relates the off-shell T-matrix to the T- 
matrix at an appropriately chosen on-shell kinematical 
point. Based on the low-energy theorem one can show 
that effects due to the off-shell structure of the T-matrix 
indeed can be ignored to a large extent. A necessary 
condition, that the full matrix element for the process 
is gauge invariant, is ensured by adding contact terms 
which are regular in the limit of vanishing photon mo- 
mentum. As a result one obtains rather accurate predic- 
tions from such a model in spite of its simplicity. Due to 
the fact that the SPM's satisfy the low-energy theorem, 
predictions are accurate as long as the nucleon-nucleon 
scattering amplitude varies little over an energy range of 
the order of the photon energy. In the past several SPM's 
have been developed for pp bremsstrahlung. The earli- 
est one is due to Low and Nyman Jl8|, |l9) and is based 
on a kind of power series exp ansion for the amplitude. 
This particular SPM [l9| will be referred to as 
Low-SPM hereafter. More recently SPM's were devel- 
oped by Liou, Lin and Gibson pi] , based on the explicit 
evaluation of the tree-level diagrams. The differences be- 
tween the different versions lie in the particular choice 
of the on-shell kinematics at which the T-matrix is eval- 
uated. Of particular interest for the discussion in this 
section is the SPM where the t and u Mandelstam vari- 
ables are selected to define the on-shell kinematics for 
the T-matrix J2l], ^2[ |2j| , which will be referred to as the 
tu-SPM hereafter. 

The novel one- and two-photon bremsstrahlung SPM 
amplitudes developed in the following are based on a 



first-order power-series expansion of the T-matrix around 
an appropriately chosen kinematical point (inspired by 
the Low-SPM) which is used in the tu-SPM. This hy- 
brid formulation (called here pse-SPM) can be shown to 
correspond to the Low-SPM with the exception of some 
additional terms proportional to the magnetic moment 
multiplied by derivatives of the T-matrix. On the other 
hand it would correspond to the tu-SPM if the T-matrix 
would depend linearly on the kinematical variables. The 
important advantage of our formulation is that the sum 
of the diagrams corresponding to radiation off external 
legs is already gauge invariant and one therefore does not 
have to introduce contact terms. 



The single-photon bremsstrahlung matrix 
element 



The antisymmetrized on-mass-shell T-matrix for 
proton-proton scattering can be decomposed in Lorentz 
scalars M as 



T(t,u) 



where covariants f2, are taken from the set 



= {1,7am0>, 75,757/J • 



(1) 



(2) 



Since only two of the three Mandelstam variables are 
independent for on-shell kinematics (s + 1 + u — 4m 2 ) we 
indicate only the dependence on (t, u) of the T-matrix 
and the invariant coefficients Cj(t,u). In order to arrive 
at the particular SPM (to be referred to as "pse-SPM" ) , 
which we will later extend to two-photon bremsstrahlung, 
it is essential to make a power-series expansion of the 
coefficients Cj of the T-matrix around a point (t 1 u) which 
corresponds to some average kinematics, 



Cf e {t,u) = C j (i,u) + (t-t) 



dCj{t,u) 



Of 



+ (u- u) 



dCj(t,u) 



du 



(3) 



where derivatives are evaluated on-shell pq j. In order 
to guarantee antisymmetry of the bremsstrahlung ampli- 
tude under the interchange of identical particles |^5| , the 
point (i, u) is defined according to 
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Ul2> 



where k 2 is the invariant mass of the emitted particle 
(k 2 = presently for a single real photon) . The Mandel- 
stam variables are defined as Si = q 2 . and similarly for 
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the others, where the 4-momenta 
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are given explicitly in terms of the momenta of the incom- 
ing and outgoing protons for bremsstrahlung (see also 
Fig- !))• 




FIG. 1: Feynman diagrams included in the calculation of the 
single-photon bremsstrahlung amplitude. The dotted lines 
depict the T-matrix, wavy lines - the photon and solid lines - 
the proton. 

Following ref. |2^] the pole contribution to the ampli- 
tude is constructed by adding the contributions from the 
four Feynman diagrams corresponding to emission from 
each of the external legs (see Fig. (|j])), 

M" = QAi(pi)«As(^) [Ti S^(pi ~ k) 

+ (pi + k) T 2 + T 3 S (2) (p 2 - k) 

+ r^ 2 > S^( P ' 2 + k) T 4 ] u X2 { P2 )u Xl { Pl ) (6) 

where the index on the T-matrix defines the kinematics 
at which it is evaluated. For the present SPM (adopted 
from the tu-SPM) the T-matrix is evaluated at an on- 
shell point defined by the same values for the (t, u) vari- 
ables as arc appropriate for the off-shell T-matrix and 
can be read from the Feynman diagrams Fig. (Q). Ex- 
pressed in terms of the momenta of the in- and out-going 
protons these are, 



Ti = T{u V2 ,t 2 ) 
T 3 = T(uu',h) 



T 2 = T(u 12 ,,t 2 ) 
T 4 = T(u V2 ,ti) 



(7) 



For the coefficients Cj{t,u) the power-series expansion 
Eq. (^) is used. Note that evaluation of the on-shell T- 
matrix at (u, t) implies for the energy s = 4m 2 , — i — 
u. This value is different from the value which could be 
inferred from diagrams in Fig. (fil). The usual expressions 
for the nucleon propagator, S(p)= i($ — irip + iO)^ 1 , and 
photon vertex (photon momentum k directed out from 
the vertex), 



pMi) 



-le 



2m„ 



(8) 



have been used, where j = 1, 2 denotes the particle num- 
ber and k is the proton anomalous magnetic moment. 

The amplitude given in Eq. (|^) has the correct pole 
structure by construction and, as can be easily checked, 



is gauge invariant without having to add contact terms. 
Thus the amplitude in Eq. (^|) obeys the low-energy the- 
orem Jl8| and qualifies as a SPM amplitude. 

Comparing the pse-SPM and Low-SPM (the version 
of ref. p5|]) in some more detail one finds that most of 
the terms are identical, with the exception of additional 
terms in pse-SPM of the type 



T^k v 



2k ■ p[ 



(t 2 -t)- + (u 12 ,~u)- 



dT 



u) 



dT 
du 



2k -pi 



+ (1 2) . 



(9) 



The obvious notation is used where ^ implies the terms 
in the Taylor-series expansion of the T-matrix that con- 
tain the derivatives of the coefficients with respect to t. 
It can be shown that the terms in Eq. (^|) are of order 
k and therefore the difference is beyond the low-energy 
theorem as one should have expected. 

It is also important to mention that the absence of the 
contact (internal) contributions in the present SPM is a 
consequence of the choice u, t as independent variables 
in the T-matrix. Should s, t (or s,u) be chosen instead, 
additional contact terms would be required to restore the 
gauge invariance (compare e.g., with the original SPM's 
of refs. |1 0). 



1. Results for single-photon bremsstrahlung 

In Fig. (||) the results for each of the three soft-photon 
models are compared with cross-section data obtained in 
a recent high precision experiment at KVI at 190 MeV 
incident energy [g3| . The predictions of the present pse- 
SPM appear to lie right in between those of Low-SPM 
and tu-SPM. As such the present SPM appears to be in 
rather good agreement with the data even though the 
photon energy is relatively large (about 80 MeV). 

Also shown is the analyzing power A y for scattering of 
the polarized protons. These results suggest that Low- 
SPM gives a better description of the data J23| at 190 
MeV than pse-SPM and tu-SPM. The latter models give 
close results. 



B. The SPM for two-photon bremsstrahlung 

The pse-SPM developed in the above can readily be 
extended to the case of two-photon bremsstrahlung since 
no explicit contact terms (which appear in the Low-SPM 
and in the tu-SPM) were introduced. One therefore does 
not have to deal with the complication of adding a two- 
photon contact term or discuss the modification of the 
single-photon contact term due to the presence of the 
second photon. 

To obtain the two-photon equivalent of pse-SPM (pse- 
2SPM) one proceeds in a similar manner as discussed in 
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FIG. 2: The differential cross section d 5 a/dQ,idQ2d8-y (upper panel) and the analyzing power A v (lower panel) for single- 
photon bremsstrahlung as functions of the photon angle. Prediction of different SPM's (see text) are compared with recent 
data obtained at KVI [E3| at a beam energy of 190 MeV. The angles of the outgoing protons are kept fixed at 8\ = 8°, 82 = 16° 
(left panel) and 81 = 12°, 82 = 18° (right panel) for coplanar geometry. 




FIG. 3: Feynman diagrams included in the calculation of the 
two-photon bremsstrahlung amplitude. 



the previous section. The T-matrix is written in terms 
of a power-series expansion around the point of average 
kinematics as given in Eq. (||), where k 2 now equals m 2 ^, 
the invariant mass squared of the two-photon system. 
This particular point is used to preserve antisymmetry of 
the matrix element (see ref 
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The amplitude can be 
constructed by adding the contributions of all diagrams 
where the two photons (with momenta q\ and (72, k — 
Qi +92) are attached to the external legs in all possible 



permutations (see Fig. (||)), 

M{seqy =u K {p' 1 )uy 2 {jp' 2 ) 
Ti S {1) {pi - k) T^ 1 ) S {1) {pi - qi) 
+all possible permutations u\ 2 (jp2)u\ 1 (pi) , (10) 

where only the first diagram of Fig. (^) has been writ- 
ten explicitly. The (t, u) variables specifying the on- 
shell point at which the T-matrix is evaluated can be 
easily expressed in terms of the external momenta for 
each diagram. However, instead of the true T-matrix the 
power-series expansion Eq. (^|) is used. It can be ver- 
ified that this amplitude satisfies gauge invariance, i.e. 
qifj,M(seq)^ u — q2 V M{seqY u = 0, for the case of radia- 
tion off the two-proton system. The amplitude therefore 
obeys the low-energy theorem for the two-photon emis- 
sion |||. 

Results for the two-photon bremsstrahlung will be pre- 
sented in a later section together with the results for the 
virtual-pion amplitude. 

III. PION PRODUCTION 

The importance of short-range physics for ir° produc- 
tion in pp scattering was addressed in many references, 



e.g. I?, 10, 1^, |27j. It was shown that this process is very 
sensitive to the short-range component of the NN interac- 
tion which in turns reflects in very strong off-shell effects 
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in the T- matrix describing the pp rescattering in the 
final state. This partial wave gives the most important 
contribution near the pion-production threshold. As a 
result of the strong off-shell effects the so-called direct 
pion production is suppressed, and other contributions, 
such as the nN rescattering, heavy-meson exchanges, etc. 
become crucial to obtain agreement with experiment. 

In this paper we have opted for simpler, more 
phenomenological approach which still relies on the 
same on-mass-shell T-matrix as used in the photon- 
bremsstrahlung calculations. The virtual-pion emis- 
sion and the sequentional two-photon bremsstrahlung 
are thus described in equivalent models. The pion- 
production amplitude is calculated using radiation off 
external legs only while evaluating the T-matrix at a suit- 
ably chosen point corresponding to on-shell kinematics. 
We use a general pion-nucleon vertex 



7 \X + 



2 m r 



(11) 



where \ specifies the admixture of pseudo-scalar (PS) 
coupling and k is the pion (outgoing) momentum. The 
PS component in the vertex is included to effectively ac- 
count for the reaction mechanisms which are not explic- 
itly present in the model. This issue will be elaborated 
further on in this section. 

As argued in || the energy dependence of real pion 
production indicates that the final-state interaction be- 
tween the emerging nucleons should be accounted for cor- 
rectly. For this reason the T-matrix is evaluated at an 
on-shell point corresponding to the same energy, s, and 
the same ratio R — t/u as is appropriate for each of the 
four diagrams. In should be noticed that for on-shell 
kinematics the ratio R is directly related to the proton- 
proton scattering angle. The amplitude which can be 
read off the diagrams in Fig. (|l|) (where the photon line 
is replaced by the pion one) has the form 



M« = T 1 sW(p 1 -k) rW{k) + 
r«(fc) S ( - 1 \p' l +k)T 2 + T s S^{p 2 
+ Yf\k) S^(p' 2 + k)T 4 



k) ri 2 )(fc) 



(12) 



with 



Ti = T(s f ,Ri = t 2 /u V2 ) , T 2 = T(s l ,R 2 = t 2 /u 12 >) , 
T 3 = T(s f ,R 3 = h/u 12 >) , T 4 = T(si, i? 4 = ti/u V2 )(13) 

The nucleon spinors in Eq. (|lj) are omitted for brevity. 

In the process of calculating the pion-production cross 
section we noticed that special care should be paid to the 
representation of the T-matrix. Initially the calculation 
was performed by expanding the T-matrix in the usual 
set of five Lorentz covariants given in Eq. (Q). 

Changing the ratio \ of the PS and pseudo- vector (PV) 
couplings by a mere 0.1 would change the real-pion pro- 
duction cross section by about one order of magnitude. 
This extreme and unrealistic sensitivity could be traced 



back to the unrealistically large coupling to negative- 
energy states in the pp system at small energies which 
is introduced with this particular choice of covariants. 
To avoid the aforementioned problem we have there- 
fore introduced another set of Lorentz tensors, chosen 
such that, when sandwiched between large components 
of positive-energy spinors, they reduce to the five oper- 
ators usually taken in a non-relativistic formulation (see 
for example RU B3, B3]), 



(!) . ^( 2 ) 



,*(<?« +a^)-n,S 12 (t),S 12 (u) , 



?(2) 



(14) 



where n = (pi —p 2 ) x (p'i—p' 2 ), t = q t , u = q u , the tensor 
is given by S\ 2 (p) = 3 crW -p -p — u^- 1 ' ■ a^ 2 ', and the 
hat denotes a unit vector. Emphasizing the dependence 
on s and R of the Lorentz-invariant coefficients Cj the on- 
shell T-matrix is expressed in a similar way as in Eq. (jl]) , 



T(s,R) = J2Cj(s,R)n j W -WW . 
A possible choice for the covariants f2j is 

^4 = 75$fc , ft 5 = ^2 + ^1 , 

defined in terms of the orthogonal four-vectors, 



(15) 



(16) 



Q? = (p?+p£)/w, 
Q£=JW-(gt-Q.Ktf), 

K(€ - (in ■ Qs)Q» 4 

^ ap (QsMQk)*(Q P ) f 



Qp 

Qn 



(«u • Qk)Qi) 



(17) 



1, 



normalized such that Q 2 S = —Q\ = —Q 2 p = —Q'. 
where W 2 = s = (pi + p 2 ) 2 and e^ uap is the fully an- 
tisymmetric Levi-Civita tensor. The momenta qt and q u 
are chosen according to Eq. (j^). It is straightforward to 
show that in the pp-CM system the matrix elements of 

flj ■ fl^ 2 ) for the large components of the spinors are 
indeed a linearly independent combination of the non- 
relativistic operators given in Eq. ( |I4| ) . The fifth term in 
Eq. (16), for example, is the only one that, in the non- 
relativistic reduction, contributes a term like the third 
one in Eq. (14). In addition the matrix elements be- 
tween large and small components vanish. This set of 
five operators is not unique; any linear combination of qt 
and q u could have been used to define Qfc an d further- 
more ill = 1 is also a valid choice. We have checked that 
any of these ambiguities have only minor effects on the 
calculated pion-production amplitude. 

With the covariants in Eq. (|l6|) the cross section for 
real-pion production still depends on the ratio x but in 
a much more gentle way. By varying x an d keeping a 
realistic nN coupling constant, the experimentally mea- 
sured cross section can be reproduced with x = 1-05. For 
X = 0, corresponding to the pure PV coupling, the cross 
section is about a factor 10 larger than experiment and 
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appears to be independent of the choice of covariants, 
Eq. (||) or Eq. (|l6|). The latter feature is probably due 
to the fact that for a PV coupling the contribution of 
negative-energy states is suppressed. 

To understand the sensitivity of the cross section to 
the PS component in the ttN vertex we can rewrite the 
pion-production amplitude in the form 



M n = M[ PV) + Mi cont) . 
where the purely PV contribution is 



(18) 



2m,, 



1 - 



k-pi — rn^/2 



75 



1 



k-p[ +ml/2 



1 



m 



k ■ p 2 — ml/2 



(2) rp 



k-p' 2 +ml/2 
and A4^° ntS> has a form of a 5-point contact vertex, 



(19) 



M 



{ - t] = -<T^)^r( Tl ^ 1) +7f ) 

V 1 + x' 2m v V 



~ T 3 75 



(2) 



(2) rri 

r 5 t 4 



(20) 



The important observation is that in the soft-pion 
limit (fc — > 0) the amplitude fulfills the re- 

quirement of the chiral symmetry |j5| (to be precise, 
the limit in Eq. (fi13) shoud be taken in the order: 
lim m7r ^o{hmg^ ...}). In this case, of course, T\ = T 2 = 
T3 = T4 is the T-matrix calculated in kinematics of on- 
shell pp scattering. However the pion mass is finite, 
and the T-matrices in the dominant diagrams, T\ and 
T3, corresponding to the pion emission off the initial 
legs, in fact should enter far off-shell (for example, the 

corresponding off-cnergy-shell CM momentum is about 
1/2 

(m„-m p ) ~ 370 MeV). Due to a strong off-shell depen- 
dence of the T-matrix (see, e.g. [^7[ |33|) a sizable re- 
duction of the cross section calculated only with Ai^ V ^ 
should occur. The contact term in Eq. ( p0| ) effectively 
accounts for this effect, as well as the other important 
mechanisms which act in the opposite direction, such as 
(off-shell) ttN rescattering and heavy-meson exchanges, 
and possible genuine contact vertices in the underlying 
Lagrangian. We will consider \ as a phenomenological 
parameter of the model. 

Our approach has a certain similarity to the soft-pion 
model of ref. |36] , where the authors applied yet another 
method to account for off-shell effects, and found a rea- 
sonable agreement with the data available at the time. 

The angle integrated cross section for real-7r° produc- 
tion is plotted in Fig. (^) as function of the relative ki- 
netic energy in the final two-proton system, defined as 
T12 = M12 — 2m„, where M12 is the invariant mass of the 
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FIG. 4: The angle-integrated cross section for real-pion pro- 
duction is plotted versus the relative kinetic energy T12 of 
the outgoing two-proton system. Results are shown for two 
incoming proton energies and for a selection of 7r-N relative 
kinetic energies. 




FIG. 5: The integrated cross section for real-pion production 
as function of 77 which is related to the beam energy through 
Eq. (|c|). The data are from §|. 



two-proton system. It can be seen that to a large extent 
the cross section is independent of T23 = A/23 — m v — 
and falls off roughly proportional to 1 /T12 in accordance 
with @. 

The total 7T° production cross section as function of ?/ 
(see Eq. ( |C5|) for the definition) is compared to the data 
of H in Fig. (||). It is seen that the calculation agrees well 
with the data in both magnitude and energy dependence. 
As such we conclude that this simple model is able to give 
a reasonable estimate of the pion-production cross section 
and will thus use it also in the calculation of virtual-pion 
production discussed in the following section. 



A. Virtual pions 

The amplitude for two-photon emission mediated by 
a virtual pion can be factorized in two terms. The first 
is the amplitude for virtual-pion production, Ai-^, which 
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is identical in structure to the one for real pions. The 
second term describes the decay of the virtual pion. The 
amplitude in question now reads 

■ 2 

J^/l^UA = %e £z[T7 e^ al3 qi q 2 

(k 2 — m 2 + i0)m7r 

xM 1T {p[,p' 2 ;pi,p2) , (21) 

where k = q± + q 2 is the momentum of the virtual pion 
and (7tt 77 ~ 0.0375 is the tt° — > 77 decay constant. The 
total amplitude for the pp — > ppr/"/ process is 

M n.v = M »»( seq ) + M^(n) , (22) 
the sum of the amplitudes given in Eq. @ and Eq. (§!]) 

IV. RESULTS 

The calculations presented in this section are done for 
an incoming proton energy of 280 MeV in the Lab system 
which is just below the pion-production threshold (2m,r + 
ml/2m p ). 

In Fig. the cross section for two-photon production 
is plotted for certain exclusive kinematics. We have opted 
to use the Dalitz coordinates (see Appendix B for a more 
detailed discussion) for expressing the differential cross 
sections, specified by i) T 12l the relative energy between 
the two protons; ii) T13, the relative energy between a 
proton and the sum-momentum of the two photons (equal 
to the momentum of the virtual pion); iii) the Euler an- 
gles of the plane spanned by the two protons and the vir- 
tual pion with respect to the incoming beam direction, 
i.e. (0k, 4>k), and where the third angle is trivial due to az- 
imuthal symmetry. This is supplemented by the angles 
(0 77 ,0 77 ) specifying the orientation of the two-photon 
relative momentum in their cm. frame and the invari- 
ant mass to 77 of the two-photon system or -equivalcntly- 
the virtual pion. These coordinates can be used for the 
sequential two-photon emission as well as for the virtual- 
pion process. We have used these coordinates instead 
of the traditionally adopted ones in bremsstrahlung for 
a few reasons. Firstly, the phase space factor is a very 
smooth (in many cases independent) function of the kine- 
matical variables and differential cross sections thus di- 
rectly reflect the magnitude of the underlying matrix ele- 
ment. Secondly, the absence of divergencies of the phase 
space factor allows for a straightforward evaluation of 
(partially) integrated cross sections. Thirdly, these coor- 
dinates uniquely determine the kinematics of the event 
while in polar coordinates a kinematical solution is not 
always uniquely defined (this happens in very selected 
parts of phase space only). It should be noted that T23 
is related to T12 and T13 by a simple algebraic relation. 

In the figures the two-photon cross sections due to 
the intermediate virtual-7r° mechanism and the sequen- 
tial two-photon emission are indicated separately. There 
is a strong interference between the two contributions, 
the total (from adding the amplitudes, labelled 'tot+' 



in the figures) is larger than the sum of the individual 
cross sections. The importance of interference for the to- 
tal cross section implies that the cross section is sensitive 
to the relative phase between the amplitudes of the un- 
corrected and the virtual-pion two-photon emission pro- 
cesses. To show this, we also plot the cross section for the 
case in which the virtual-pion matrix element has been 
arbitrarily, only for display purposes, multiplied with a 
minus sign (i.e. changing the relative sign in Eq. (|22] ) 
eventhough the sign given there is correct). This calcula- 
tion, labeled 'tot-' in Fig. @, gives rise to a much smaller 
total cross section. It should be noted that changing the 
sign does not affect the real-pion production cross sec- 
tion since it is independent of the sign. From Fig. (|^) 
it can be seen that the angular distributions depend on 
the phase of the virtual-pion contribution, however the 
largest effect of changing the sign shows up in the overall 
magnitude. We have checked that this is the case in a 
large region of phase space and the distributions shown 
here can be regarded as typical. 

Another aspect which can be seen from Fig. (||) is that 
the angular distributions of the sequentional two-photon 
emission process show pronounced structures. This is to 
be expected as the single-photon bremsstrahlung angular 
distribution shows pronounced peaks which are due to 
the quadrupole nature of the electric radiation and the 
interference with magnetic radiation. The virtual-pion 
mechanism has a rather featureless distributions, due to 
the fact that the two photons couple to the quantum 
numbers of the pion, J 77 = 0~ . 

It is also apparent from Fig. (|J) that -unfortunately- 
one cannot point to a particular feature in the angular 
distribution which is especially sensitive to the virtual- 
pion contribution. There are no quantum numbers that 
distinguish this process from the bremsstrahlung contri- 
bution. 

For the above reason, and also because cross sections 
are -in general- small for two-photon emission, we inves- 
tigate whether the virtual-pion signal can also be seen in 
less exclusive kinematics where certain angles have been 
integrated. Since, as remarked before, the virtual pion 
contribution seems to give primarily rise to an overall in- 
crease of the cross sections we have performed a simple 
integration of the differential cross section. 

The squared matrix element for pion emission is in- 
versely proportional to the relative energy in the final 
pp system |||. One thus expects that the virtual-pion 
process is most pronounced for the lowest values of T\ 2 . 
This is indeed supported by our calculations as shown in 
Fig. (0) , where the difference between the full calculation 
(labeled 'tot+') and the sequential two-photon process 
strongly depends on T12 and hardly on T23. In Fig. (0) 
all angles have been integrated. 

The unambiguous signature of the virtual-pion contri- 
bution is that it increases the closer one approaches the 
real-pion pole. This can clearly be seen from Fig. (j|) 
where the cross section is shown as function of m 77 at 
fixed T12. All other variables, i.e. all angles and also T23, 
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E=280MeV; M =110MeV; T 13 =20MeV; T 23 =24MeV 
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FIG. 6: The dependence of the differential cross section (fia/dmjjdM^dM^dQjjdQkdfj)! on the two-photon angles. The two- 
photon bremsstrahlung result is given by the dense-dotted curve, the pure virtual-pion process by the sparse-dotted curve, the 
full result by the drawn curve, while the dashed curve gives the coherent sum when changing ad-hoc the sign of the virtual-pion 
contribution. 
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FIG. 7: The cross section for two-photon production, inte- 
grated over angles, is plotted versus the relative kinetic en- 
ergy T12 of the final two-proton system. Results are shown 
for a selection of 7r-N relative kinetic energies. The meaning 
of the curves is the same as in Fig. (|^) . The end points of the 
curves are determined by kinematics. 



are integrated. Even this rather inclusive cross section 
shows a clear sensitivity to the interference between the 
sequential and the virtual-pion two-photon emission pro- 
cesses. 
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FIG. 8: The cross section for two-photon production, in- 
tegrated over all variables at fixed relative kinetic energy 
T12 = 3 MeV of the final two-proton system, is plotted as 
function of the two-photon invariant mass. The meaning of 
the curves is the same as in Fig. ([]). 



V. SUMMARY AND CONCLUSIONS 

In this work we have shown that the two-photon 
bremsstrahlung offers the interesting possibility to 'mea- 
sure' subthreshold pion production in pp scattering. It 
allows for studying pion production in kinematics which 
is not accessible in the pp —> ppir reaction. In addition, 
the phase of the virtual-pion process with respect to that 
of sequential two-photon emission can be investigated. 

To account for the sequential two-photon emission pro- 
cess, which is an important background, a novel soft- 
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photon model (called pse-SPM) is developed. This model 
is tested in a calculation of single- photon bremsstrahlung, 
and is shown to give accurate results for cross sections. 

Calculated exclusive cross sections of the pp — > ppjj 
reaction are in general small, however sensitivity of the 
cross sections to the virtual-pion signal remains even for 
rather inclusive cross sections. 
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APPENDIX A: KINEMATICS FOR 
TWO-PHOTON PRODUCTION 

For the reaction N + N — > N + N + 7 + 7 the momenta 
are denoted by pi,p 2 , p\ , p' 2 , Qi , 92 (see Fig. (§) ) . Energy- 
momentum conservation reads P1+P2 = Pi+p' 2 + qi + q2- 
The cross section is 

da = -S- j \A\ 2 {2vf5\p 1 +p 2 -p' 1 -p l 2 -q 1 -q 2 ) 

d 3 p[ d 3 p' 2 d?qi d 3 q 2 
X (2ir) 3 E[ {2tt) 3 E' 2 (27r) 3 2ei (2^) 3 2e 2 ' 

where A = -M^e^e^ is the invariant amplitude, e\ and 
e 2 are the polarization vectors of the photons, £\ = \qi\, 



£ 2 = \q 2 \, and j = J (p x ■ p 2 ) 2 



,4 - 



m p \p La b\ in the 

laboratory frame where p 2 — (m p ,0). Using the identity 
J & A {qi + q 2 — k)d A k = 1 the cross section is put in the 
form 

771 4 f 

do = t^t J \A\ 2 5 i {p 1+ p 2 -p' 1 -p'. 2 -k)I 11 



d 3 p[ d 3 p' 2 4 



(Al) 



To calculate this integral in an arbitrary frame we intro- 
duce the relative and total 4-momenta of the photons 



k = q 1 + q 2 , 



1 = \{<li - 92), 



qi = -k + l, q 2 = ~k-l. 



(A2) 



The Jacobian of the transformation from q\, q 2 to I, k 
is unity, and after removing the trivial <5-function we get 

fu , % dH 

77 = J ,(6 l w+B * w " ko) S^^r 
\i\ 2 



4(|/]feo — \k\h cos# 77 ) 



-dfl 



77 ) 



(A3) 



with dQy-y = sin# 77 G?# 77 d0 77 , where we introduced the 
polar and azimuthal angles # 77 and 77 between the 3- 
vectors k and I. For real photons (qf = q 2 = 0) one can 



show that k ■ I = and 4 I 2 - 



■mz 



0, where m 2 = fc 2 is 



77 YViiV^i^ .IV^^ 

the invariant mass of the two-photon system. Expressing 
now Iq in terms of the 3-momentum \l\ we obtain 



2a/1 



cos^ 



77 



k 2 sin 2 



with kn — 



Eq. ( |A3| ) can be simplified to 

|Z| 3 



77 + k 2 . The two-photon phase space 
iplified to 

1 n 3 

(A4) 



k m 2 



As a last step the integration over ko in Eq. (Al) is 
replaced by an integration over the two-photon invariant 
mass using kgdko — m 77 <im 77 . We obtain 



da = 



2ml 



I A j 2 J (m 77 ) Ijj m 77 dm 



<"rr > 



(A5) 



(27T)8j 

where we introduced the 3-particle phase-space integral 



J(m 77 ) 



S 4 (pi 



■P2 



1 

Pi 



P 2 ~k) 



d 3 p[ d 3 P ' 2 d 3 k 
' E[ E' 2 2k 



(A6) 



In the one-photon bremsstrahlung the similar integral is 
traditionally evaluated in polar coordinates (see, for ex- 
ample, Eq]) leading to the cross section of the type shown 
in Fig. (g). For the two-photon bremsstrahlung in the 
present paper we will use the Dalitz coordinates instead, 
as discussed in Appendix B. 



where J 77 is the two-photon phase-space integral defined 
as 



J 77 = J 5 A {qi +q 2 - k )-^ 



d 3 qi d 3 q 2 



2e 2 



APPENDIX B: DALITZ COORDINATES 

To evaluate the phase-space integral in Eq. ( |A6| ) we 
choose the CM frame where pi + p 2 = 0, and carry out 
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the integration over p 2 . Introducing s — (pi + P2) 2 we 
obtain 



J(m 7 



where E 2 = m 2 + p 2 — m 2 + (p[ + k) 2 = m 2 + pf 
k 2 + 2p[k cos 0i3. Using 



d 3 p[ d 3 k = pi E[ dE[ dtti k k dk dn k 

= p[ E[ dE[ d cos dcj)\ k k dk d cos Ok dtfrk 

and integrating over cos #13 using the (5-function we ob- 
tain 

J(m 77 ) = —d cos Ok dcf>k d(f>\ dko dE[ . 

Defining invariant masses 

M 2 2 = (T 12 + 2m p ) 2 = (pi+p^) 2 

= s — 2y/sk + »™ 77 , 

M 2 3 = (T 23 + m p + m^f = (p' 2 + k) 2 

= s- 2y/sE[ +m 2 p , 

we can cast the integral J(m 77 ) in the form 
J(m 77 ) = -^-dcosQk d(f>k depi dM\ 2 dM 2Z . 



(Bl) 



(B2) 



Here the angles 4>k, Ok, 4>i describe in the CM frame the 
orientation of the plane in which the momenta lie of the 
outgoing two protons and the two-photon system (the 
virtual pion) with respect to the incoming beam. Specif- 
ically, the angle 4>\ is defined as the azimuthal angle of 
the momentum p[ in the frame, where momentum k is 
along OZ axis and the OX axis lies in the plane spanned 
by the beam and fc. The angle Ok is taken as the angle 
in the CM frame between the incoming momentum and 
k while <pk is a trivial azimuthal angle. In the CM frame 
the momenta can now be expressed as 

k = (k x , k y , k z ) = k(sinO k , 0, cos 9k), 

Pi = [Plx, Ply, Plz) = 

= p[ (cos Ok sin 0i3 cos 4>\ + sin Ok cos #13, 
sin 0i3 sin 0i , — sin Ok sin 0i 3 cos (f>i + cos Ok cos 013), 

P2 = ~P{ - k - 

The magnitudes of fc, p[ and p 2 are determined through 
the energies of the two-photon system and nucleons 



fc 



E> 



'12 



,2 



2yfs 



s/s-E' 1 -k , 



and the angle between p[ and k can be expressed as, 
cos 0i3 = {E 12 - E' 2 -k 2 + m 77 )/2fcp'i- 

So far the momenta were defined in the CM frame. The 
boost to the Lab system is specified by the velocity V = 
Pi I {fnp + Ei) and the Lorentz-factor 7 = (m p + Ex)/ y/s. 
The Z components of the vectors in the lab can now be 
expressed as 



k Lab 

/Lab 
Plz 

lLab 
P2z 



j(k z + Vk ), 
jip'^ + VEi), 

pl~k L z ab -p[ L z ab , 



while the X and Y components do not change. 



APPENDIX C: TOTAL CROSS SECTION FOR 
PION PRODUCTION 

Using the Dalitz coordinates, the cross section for real 
pion production is written as 



da 



(27r)5j 



\M„\ 2 J(m„) 



(CI) 



where particle 3 is associated with the pion and J(m w ) 
is given in Eq. ( |B2[ ) where m 77 — » m n . 

Let us first make a tentative assumption that the am- 
plitude is a constant and integrate the cross section in 
Eq. JCl]) over angles, 



da = Wjsfs lM * |2 dM * 2 dM ™ ■ (C2) 

For the total cross section we integrate over invariant 
masses 



A 2 

m,,7r 
P \MA 2 I{s), 



{2nfjs 

where we introduced the integral 



(C3) 







^ = A 









[f 2 



= 2 



M 
M 



yJ(M 2 2 - 4m 2 ) 



12 



xyj(s-ml- M 2 2 ) 2 - 4mlM 2 2 dM 12 ,(C4) 

with the lower and upper limits M^ m — 2m p and 
j[frnax _ ^ _ m ^ respectively. At the upper 

limit the pion 3-momentum in the CM system, k 2 — 
(M"i 2 2 - (y/s - m^) 2 ) (M\ 2 - (y/s + inn) 2 ) /4s , vanishes 
while it reaches a maximum at the lower limit, 

k 2 max = ml V 2 = b(b + im^)/4s (C5) 
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which defines the conventionally used variable rj ||] in 
terms of b — (y/s — m^) 2 — Am 2 Introducing the rela- 



tive pp 3- momentum Ap 2 — (Mf 2 — Am 2 ) the phase space 
integral Eq. can also be put in the familiar form 

I(s)=J16\p\\k\V^W2df- 

With the substitutions t = (M 2 2 - Am 2 )/b, a = Am 2 

and c = (y/s + m^) 2 — Am 2 the phase space integral 



Eq. (C4) can be cast in the form 



1(a) 




t(t- \){t-c/b) 
(t + a/6) 



(C6) 



This shows that the total cross section is roughly propor- 
tional to b 2 and thus proportional to rf at low energies, 
in contrast to the data. It has been argued in || that the 
assumption made that the matrix element is constant is 
not valid. The pion-emission process is strongly affected 
by the final-state NN interaction at low energies, the ef- 
fect of which is roughly proportional to 1/p 2 = A/tb. 
Including this factor in the integrand gives a total cross 
section proportional to rj 2 , in rough agreement with the 
data (above rj = 0.2). The approach that we chose for 
evaluating the amplitude and comparison with experi- 
ment arc described in detail in sect. III. 
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